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ABSTRACT 

The net primary productivity (C) or the 
annual rate of carbon accumulation per 
unit ground area by terrestrial plant 
communities is the difference of gross 
photosynthesis (A g ) and respiration (R) 
per unit ground area. Available field 
observations show that R is a large and 
variable fraction of A g , although it is 
generally recognized that there are 
considerable difficulties in determining 
these fluxes, and thus pose challenge in 
assessing the accuracy. Further 
uncertainties arise in extrapolating field 
measurements (which are acquired over 
a hectare or so area) to regional scale. 
Here, an approach is presented for 
determining these fluxes using satellite 
and ancillary data to be representative of 
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regional scale and allow assessment of 
interannual variation. A g has been 
expressed as the product of radiation use 
efficiency for gross photosynthesis by an 
unstressed canopy and intercepted 
photosynthetically active radiation, 
which is then adjusted for stresses due to 
soil water shortage and temperature 
away from optimum. R has been 
calculated as the sum of growth and 
maintenance components (respectively, 
R g and R m ).The R m has been determined 
from nitrogen content of plant tissue per 
unit ground area, while R g has been 
obtained as a fraction of the difference 
of A g and R m . Results for five 
consecutive years (1986-1990) are 
presented for the Amazon-Tocontins, 
Mississippi, and Ob River basins. 

INTRODUCTION 

The net primary productivity or the rate 
of carbon accumulation by terrestrial 
plant communities per unit area (C) is 
the difference of gross photosynthesis by 
the canopy (A g ) and autotrophic 
respiration of the stand (R): 

C = A g - R (1) 

A rather attractive approach for 
determining C was proposed by 
Monteith (1977) as: 



C = RUE * IPAR (2) 

where IPAR is the intercepted 
photosynthetically active radiation and 
RUE is the radiation use efficiency. This 
approach has been used to determine C 
at varied spatial scale (e.g., Ruimy et al., 
1994). 

Seasonal variation of A g is determined 
primarily by IPAR and the maximum 
rate of photosynthesis by leaves 
(Monteith, 1977; Choudhury, 2000a), 
while seasonal variation of R is 
determined by changes of leaf, stem and 
root biomass and their specific 
respiration coefficients, and these 
coefficients have been observed to 
increase with tissue temperature 
(Choudhury, 2001a). Because of these 
differing sensitivities of A g and R to 
biophysical factors, it is desirable to 
calculate these two fluxes separately. 

MODELING AND INPUT DATA 

Gross Photosynthesis 

The approach used to determine A g can 

be expressed as follows: 

A g = RUEg * IPAR * Ei (W) * E 2 (T) (3) 
where RUE g is radiation use efficiency 
for gross photosynthesis by an 
unstressed canopy, and Ei (W) and E 2 
(T) are prescribed functions of available 
soil moisture (W) and foliage 
temperature (T). These functions 
account for stresses due to soil water 
deficit, and deviation of foliage 
temperature from the optimum range for 
photosynthesis. 

Details about the calculation of A g and 
input data requirements can be found in 
Choudhury (2000a), and briefly 
described here. The variation of 
photosynthesis by a leaf with irradiance 
has been described by a non-rectangular 
hyperbola. The maximum rate of leaf 
photosynthesis in this non-rectangular 
hyperbola depends upon the leaf 
nitrogen content, which varies within the 
canopy. The radiative transfer equation, 
with boundary condition of direct and 


diffuse PAR incident on the canopy, is 
solved to obtain irradiance on leaves, 
and interception by the canopy. The 
single scattering albedo of the radiative 
transfer equation, which controls 
scattering of irradiance by leaves, also 
depends upon the nitrogen content of 
leaves. The sunlit leaves receive both 
direct and scattered irradiance, while 
shaded leaves receive only the scattered 
irradiance. The instantaneous canopy 
photosynthesis is obtained by 
numerically integrating photosynthesis 
by sunlit and shaded leaves weighted by 
their fractions at any depth within the 
canopy. The instantaneous intercepted 
PAR, and canopy photosynthesis are 
integrated at 1 5 min. interval throughout 
the day-light period considering diurnal 
variation of incident direct and diffuse 
PAR to obtain the daily total A g for an 
unstressed canopy and IPAR. The 
available soil moisture needed to 
account for the effect of soil water stress 
on A g has been determined by solving 
water and energy balance equations 
(Choudhury and DiGirolamo, 1998). 

The input meteorological data (e.g., 
incident PAR, air temperature, 
precipitation) have been derived from 
radiances measured by sensors on board 
polar orbiting and geo-stationary 
satellites (e.g., NOAA, GOES), obtained 
as output of four dimensional data 
analysis (specifically the frictional 
velocity, needed to solve the energy 
balance equation) and as a merged 
product of ground observations and 
satellite-based measurements 

(specifically, precipitation used in the 
water balance calculations). The 
fractional interception of incident 
photosynthetically active radiation has 
been calculated from soil adjusted 
vegetation index (SAVI), which is 
determined from visible and near- 
infrared reflectance derived from 
measurements by the advanced very 
high resolution radiometer (AVHRR) on 


board the NOAA satellites after 
correcting from sensor degradation and 
atmospheric affects on the satellite- 
observed radiances (Choudhury et al., 
1994; Choudhury and DiGirolamo, 
1998). 

Respiration 

Respiration (R) has been calculated 
according to the functional, two- 
compartment model of McCree (1974), 
as being the sum of growth and 
maintenance processes (denoted, 
respectively, as R g and R m ). The R g has 
been expressed as: 

Rg = (1"Yg) { A g - R m } (4) 

where Yq is the growth conversion 
efficiency (i.e., mole of carbon 
synthesized per unit mole assimilated). 
Chemical composition of plant tissues 
determines the magnitude of Yg and it is 
not directly affected by stress. A 
synthesis of the published values of Yq 
for different plant species give 0.69 for 
woody vegetation (shrubs and trees) and 
0.74 for herbaceous vegetation (natural 
grasses and crops), with a coefficient of 
variation of about 12%. 

R m has generally been observed to vary 
strongly with tissue temperature, and it 
has been expressed in terms of its value 
at a specified mean daily temperature 
(T 0 , °C), namely, R m (T 0 ), as: 
R m (T)=R m (T o )Q 10 {(t - To)/10) I o (AT/2) (5) 
where Qio is the temperature response 
coefficient (the value of which 
determines the magnitude of change in 
Rm (T) due to 10°C change in 
temperature), I 0 (y) is the modified 
Bessel function of argument y, and AT is 
the diurnal range of tissue temperature. 
Observations give Qio to be about 2.0 
around 20°C, while it increases for 
lower, and decreases for higher 
temperature. 

Available data for R m (T 0 ) at any 
prescribed T c show much variability 
(more than a factor of three) for different 
tissues (roots, stem, and leaves) of a 
particular species, and this is also the 


case for the same tissue for different 
species. Observations for different 
tissues and species analyzed by Ryan 
(1991) show that while R m (T 0 ) varied 
by three orders of magnitude, 58% of the 
variability of Rm can be explained as 
being due to the associated variability of 
the nitrogen content (N) of the tissue (as 
determined by linear regression 
analysis). Such a relationship between 
R m and N can be reconciled by noting 
that protein turnover in plant tissues is 
considered to a major determinant of R m . 
For a given nitrogen content, R m (T 0 ) has 
been observed to be about a factor two 
higher for roots as compared to that for 
leaves, which has been suggested to be 
due to the work done by roots to extract 
ions from soil. Observations also show 
that, for a given nitrogen content, R m 
(T 0 ) of leaves for species growing in 
boreal climate is about a factor of two 
higher than that for species growing in 
temperate climate, which has been 
suggested to be due to adaptation to 
growing in noticeably different 
temperature. 

The equation for R m (T o =20) (in units of 
mol CO 2 m' 2 d' 1 ) for species growing in 
a temperate region used in this study is: 
R m (T 0 ) = 0.21 ( N f + N s + 2 N r ) (6) 

where Nf, N s , and N r are, respectively, 
the nitrogen content per unit ground area 
of live foliage, stem, and roots (mol N 
m" 2 ). The first, second and the third term 
on the right hand side of Eq. (6) 
provides, respectively, maintenance 
respiration of foliage, stem and roots. As 
noted above, the cost of maintaining 
roots is higher than that for foliage. 

The input data required for calculating 
R, in addition to those required for A g , 
are the diurnal range of temperature 
plant tissues (AT) and L 0 . In this study, 
AT has been approximated to be that for 
air temperature, which was determined 
by combining surface and satellite data. 
L 0 has been determined from the 
AVHRR data. 



RESULTS AND DISCUSSION 
The results of calculations for 60 
consecutive months (1986-1990) are 
presented in Table 1 for the Amazon, 
Mississippi, and Ob River basins. None 
of the model parameters have been 
determined by calibration against field 
measurements of these fluxes or to 
match fluxes. Respiration is calculated to 
be about 73%, 63%, and 66% of gross 
photosynthesis of these three river 
basins, respectively. Particularly 
noticeable is the interannual variation of 
two of the three river basins. Rather low 
gross photosynthesis of 60 mol C m' 2 yr 
was calculated for the Mississippi 
River basin for the year 1988, due 
largely to drought year for most of the 
basin. Maintenance respiration is 
calculated to be about 75% of total 
respiration. Representativeness of these 
results are being evaluated. 

Table 1. Calculated Mean and Inter- 
annual Range of Carbon Balance Fluxes 
(mol C m' 2 yr* 1 ) for the Amazon, 
Mississippi, and Ob River basins for 
1986-1990. 


A g R C 

Amazon River basin: 

Mean 194 142 51 

Range 187-200 140-146 49-59 


Mississippi River basin: 

Mean 80 50 30 

Range 60-89 47-52 13-37 

Ob River basin: 


Mean 56 37 19 

Range 51-64 30-40 13-26 
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